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d 1 INTRODUCTION iMetz et al. 20081 ). Another approach to explain galaxy ro- 



ABSTRACT 

We investigate the mean velocity dispersion and the velocity dispersion profile of stel- 
lar systems in MOND, using the N-body code N-MODY, which is a particle-mesh 
based code with a numerical MOND potential solver developed by Ciotti, Londrillo 
and Nipoti (2006). We have calculated mean velocity dispersions for stellar systems 
following Plummer density distributions with masses in the range of 10^ M0 to IQ^Mq 
and which are either isolated or immersed in an external field. Our integrations repro- 
duce previous analytic estimates for stellar velocities in systems in the deep MOND 
regime (ai,ae ^ ao), where the motion of stars is either dominated by internal ac- 
celerations (fli 3> tte) or constant external accelerations (og 3> a^). In addition, we 
derive for the first time analytic formulae for the line-of-sight velocity dispersion in 
the intermediate regime (a.^ ~ Oe ~ ao). This allows for a much improved compari- 
son of MOND with observed velocity dispersions of stellar systems. We finally derive 
the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way 
halo globular clusters that have recently been proposed as a differentiator between 
Newtonian and MONDian dynamics. 

Key words: galaxies: clusters: general- galaxies: dwarf - gravitation - methods: 
analytical - methods: N-body simulations 



The flattening of rotation curves of disk galaxies at large 
radial distances, i.e. the apparently non-Newtonian mo- 



tation curves would be an alternative theory of gravity. One 

promising alternative theory is Modified Gravity (MOG) 

which has rec ently been successfully applied for dwarf satel- 
tion, IS usually explained by invoking the otherwise unde- , . a- ^ c rn r,r,r,r, h 1 j- ^ 11 1 

' , „ , J; , , , , , ,^-r^^,s TT^^nn lite gala xies HMoffat fc Toth 2007al ) and distant globular 

tected, so called Cold Dark Matter (CDM) (Bosma 19811 : , ^ Jm- n\ rn nrvo^i l\ ,^ j_ r 1 

1 .^ , . ' " : — m, . , , . , ^ TTT clusters (IMoffat fc Toth 2007bl ). One of the most famous al- 

IHubm fc Burstem 19851 ). this hypothesis has successtully , ,• ,1 • ■ ,u n j j-c j ^ ■ j 

; ; — ; f^' r 1 ternative theories is the so-called modifaed newtonian dy- 

explamed the internal dynamics of galaxy clusters, grav- . /,j^att^\ ii 1 • 1 • ^ 1 1 1 m-i 

, , , , r 1 namics (MOND) theory, which was introduced by Milgrom 

itational lensmg and the standard model of cosmology r-,nan\ a j- x TM-r~>ATT^ ii a ^ ^ i- r • 

" If 1 1 • ■ /■nT,\ • (1983). According to MOND, the flat rotation curves ot spi- 



within the framework of general relativity (GR). Despite 
the fact that the dark matter model h as been notably suc- 
cessful on large scales l|Spergel 20031 ). dark matter parti- 



ral galaxies at large distances can be explained by a modifi- 
cation of Newton's second law of acceleration below a char- 
acteristic scale of ao — 10~^°ms~^ without invoking dark 

cles has not been detected after much experimental ef- 177-; — — „ , ... , .1, 

. matter (|Bekenstein fc Mflgrom 1984 ). 



It has been shown that on galactic scales MOND 



forts and the results of high resolution N-body simula- 
tions do not seem to be compatible with observations , . , , „_ , . 

, ^. , jtt: — : — .„„^ [77 — — .^rvrJ can explain many phenomena at least as well as CDM 

on galactic scales flKlvpm et al. 1993 : Moore et al. 1993 : jr; . - — , ^r^^^^ n 1 n 1 r. 1 1 

a „ — (|Sanders fc McGaugh 20021 ). For example, Sanchez-Salcedo 

and Hernandez (2007) studied the tidal radii of dis- 
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(HB); pavel@astro.uni-bonn.de (PK); ies in MONDian dynamics. The most serious chal- 
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l|Sanders fc McGaugh 2o"o3 ). and the merging of galaxy 
clusters, where the baryonic mass is clearly separated 
from the gr avitational ma ss, as indicated by gravitational 
lensing (IClowe et al. 20061 ). Both phenomena can be ex- 
plained in MOND if some kind of hot dark matter is as- 
sumed, perhaps in the form of a massive (~ 2eV) neutrino 
l| Angus et al. 20061 1. 

MOND has rec ently been gener alized to a general- 
relativistic version (|Bekenstein 200 j ). making it possi- 
ble to test its predictions for gravitational lensing. In 
the non-relativistic version of MOND, the acceleration 
of a particle due to a mas s distribution p is given by 
l|Bekenstein fc Milgrom 1984 ): 

V- — )a) =47rGp = V-ajv, (1) 
Co 

where ajv is the Newtonian acceleration vector, a is the 
MONDian acceleration vector, a = |a| is the absolute value 
of MONDian acceleration and /i is an interpolating func- 
tion which runs smoothly from ^j.{x) = a; at a; ^ 1 to 
IJ.{x) = 1 at a; 3> 1. The standard interpolating function 
is ^ii{x) — I ^ , but Famaey & Binney (2005) suggested 

another function ii2(x) — j^, which provides a better fit 
to the rotation curve of the Milky Way. Equation ([T} can be 
transformed into V- (/i(^)a— a]v) = 0, where the expression 
in parentheses is thus a curl field, and we may write 

^( — )a = aiv 4- V X H. (2) 

ao 

The value of the curl field H depends on the boundary con- 
ditions and the mass distribution, but vanishes for some spe- 
cial symmetries. In realistic geometries, the curl field is non- 
zero and leads to difficulties for standard N-body codes. In 
other words, the non-linearity of the MOND field equation 
makes the use of the usual Newtonian N-body simulation 
codes impossible in the MOND regime. 

Many stellar systems (e.g. globular clusters) have tidal 
radii much larger than their sizes, therefore the external field 
is approximately constant over the cluster area and the mo- 
tions of stars are not influenced by tidal effects. 

In Newtonian dynamics, a stellar system evolving un- 
der the influence of a uniform external acceleration, will, in 
the frame of the system, have the same internal dynamics 
as an isolated system. In MOND, due to the non-linearity of 
Poisson's e quation, the strong equivalen ce principle (SEP) 
is violated l|Bekenstein fc Milgrom 19841 '). and consequently 
the internal properties and the morphology of a stellar sys- 
tem are affected both by the internal and external fleld. This 
so-called external fleld effect (EFE) signiflcantly affects non- 
isolated systems and can provide a strict test for MOND. 
The EFE postulation originated from observations of open 
clusters in the solar neighborhood, which do not show mass 
discrepancies e ven if the internal accelerations are below ao 
l|Milgrom 19831 ). The EFE has several consequences, for ex- 
ample it allows high ve locity stars to escape from the poten - 
tial of the Milky Way (jFamaev et al. 20071 : IWu et al 20071 ) . 
and it decreases the velocity of satellite galaxies at very large 
radii, which is in c onflict with the asymptotically flatt ening 
of rotation curves l|Gentile et al. 20071 : IWu et al. 20081 ). The 
EFE implies that rotation curves of spiral galaxies should 
fall where the internal acceleration becomes equal to the 
external acceleration. In addition, if the EFE is taken into 



account, internal properties of Gal axies such as t he TuUy- 
Fisher relation should be changed l|Wu et al. 20071 ). 

Milgrom derived the mean velocity dispersion of stellar 
systems for two special cases of internal or external fleld 
dominated systems 

analytically, assuming that the systems are everywhere 
in the deep-MOND regime {a^,ai <C ao)- If the external 
acceleration is much larger than the internal one a^, the 
system of mass M is in the quasi Newtonian regime but with 
a normalized gravitational constant larger than the standard 
Newtonian one by a fa ctor ^ , and th erefore the line-of-sight 
velocity dispersion is jMilgriam 19861) 

ae 

where (tlos,n is the Newtonian velocity dispersion. If Oe <^ 
Ui <^ ao, the cluster is isolated and the line-of-sight velocity 
dispersion is given by 

aLOS,M2 = OA71{GMao)i. (4) 

Many systems which can be used to test MOND are not 
completely internally or externally dominated, for example 
globular clusters or dwarf galaxies of the Milky Way have 
inter nal and external acceler ations which are of the same or- 
der (|Baumgardt et al. 20051 ). Since Milgrom's relations are 
valid only for systems that have either a^ ^ Oe or at <^ 
and are in the deep-MONDian regime, one has to deter- 
mine the velocity dispersions numerically for intermediate 
cases. Milgrom found that for isolated systems (internal ac- 
celeration dominated), the mass M of a system is nearly 
proportional to the forth power of the line of sight velocity 
dispersion aios and the ratio af^^/GM must be somewhere 
between |ao and ao. But how does the velocity dispersion 
change while the system transits from the Newtonian to the 
MONDian regime? In an attempt to answer this question, 
we have performed N-body simulations and present analyt- 
ical formulae for the velocity dispersion of stellar systems 
in the intermediate MOND regime. We have calculated the 
velocity dispersion for a number of isolated systems in which 
the internal accelerations a; are in the range from a^ ^ ao 
to Oi S> ao. We also give formulae for systems with differ- 
ent strengths of external flelds. It should be noted that the 
isolated systems are in equilibrium only in the Newtonian 
case, and reach a MONDian equilibrium state after collapse. 
For non-isolated systems we start from the MONDian equi- 
librium state which is created as described in section 4.3. 
These results could be useful for comparison with observa- 
tional data of several GCs and dSph galaxies that are far 
away from the host galaxy, so that the external acceleration 
due to the host galaxy is small (a^ < O.Olao) these objects 
should therefore provide straightforward possibilities to test 
MOND. Since the external fleld affects the velocity disper- 
sion by both tidal effects and EFE, and in order to see the 
pure MONDian effects, we concentrate on systems in which 
the tidal radius is much larger than the virial radius and 
therefore tidal effects are unimportant. 

This is the flrst of a series of papers that deals with 
the numerical calculations for stellar systems. In the forth- 
coming papers, the observational constraints on mass and 
velocity dispersion of Pal 14 will be studied by Hilker et al. 
(2008) and Jordi et al. (2008). 
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The paper is organized as follows: In Section 2 we in- 
troduce theoretical predictions for the velocity dispersion 
in different regimes. In Section 3, we give a brief review of 
the N-MODY code which we use for our modelling. The 
numerical results for isolated and non-isolated systems and 
comparison with observational data are discussed in Section 
4. We present our conclusion in Section 5. 



2 LINE-OF-SIGHT VELOCITY DISPERSION 
IN DIFFERENT REGIMES 

2.1 Newtonian regime 

In Newtonian gravity, the mean-square velocity, ct^, of a stel- 
lar system of mass M is given by the following equation 
(Equation (4-80a) of Binney and Tremaine (1987)): 

2 GM 



(5) 



where rg is the gravitational radius defined as (Equation 
(2-132) of Binney and Tremaine (1987)): 



\W\ 



(6) 



Here W is t he total potenti al energy. In the case of a Plum- 
mer model (jPlummer 191lh . and if we assume an isotropic 
velocity distribution, the line-of-sight velocity dispersion be- 



0'LOS,N = 0.36 



GM 

Rh 



(7) 



where Rh is the half-mass radius. If we define the half-mass- 

GM 



radius acceleration as ah = 
relation as 



'^LOS,N 

GM 



0.035ah. 



, we can re-write the above 



(8) 



2.2 MOND regime 

In the case of MOND, and in the presence of an external 
field, the total acceleration, which is the sum of the inter- 
nal ai and extern al tte acceleration, satisfies the modified 
Poisson equation (iBekenstein &: Milgrom 19841 ') . 

V.[^(^l±^)(ai + ae)] ~ 47TGP, (9) 
ao 

where fle is approximately constant, ai = V(f> is the non- 
external part of the potential and p is the density of the star 
cluster. The boundary condition is \/<j) — a^x for r ~* oo. 
Equation[9]was postulated by Milgrom (1983) to explain the 
dynamical properties of nearby open clusters in the Milky 
Way and is an outcome of the MOND phenomenology. As an 
approximation for a spherical system one can write equation 
Has: 

Me+ai\ 

aifi{ )=aiv. (10) 



ao 



Note however that Ea llQI is only an approximate and 
effective way to take into account the external field effect 
(EFE), in order to avoid solving the modified Poisson equa- 
tion with an external source term pext on the right-hand 
side. The EFE is indeed a phenomenological requirement 



of MOND, which has important consequences for non- 
isolated systems. For example, if Oe <C a* <^ ao, then 
the dynamics is in the MOND regime, and the external 
field can be neglected. When ai ^ ae ^ ao, p, tends 
to its asymptotic value p{ae/ao) ~ a^/ao {saturation of 
the p function), and the gravitational potential is thus 
Newtonian with a renormalized gravitational constant to 
{Geff = G/p{ae/ao) « Gao/ae f llMilgrom 198 6h]V 
Recently, several papers were published using 
this formulation to take into account the EFE 



Gentile et al. 2007 



Wu et al. 2007 



Wu et al. 200: 



[Angus 20081: iFamaev et al. 20071 : iKlvpin fc Pr ada 200 
For example, in order to estimate the order of magnitude 
of the EFE, Famaey et al. (2007) and Gentile et al. (2007), 
pointed out aoa{a + ae)/{ao + a + Oe) = ajv using a simple 
/ ^-function. O t her authors replaced \ai + ae\ = a^ + a^ 
(|Angus 20081 : iKlvpin fc Prada 20081 ). A more rigorous 
treatment of EFE on galactic rotation curves was made by 
Wu et al. (2007, 2008). In the work by Wu et al. (2007), for 
the mass density of the internal system, the MOND Poisson 
equation was solved as if the system was isolated, but the 
boundary condition on the last grid point was changed to 
be nonzero. 

As a first approximation, we considered that a cluster is 
in a non-inertial frame, which free-falls with a uniform 
systematic acceleration. Since the calculation of is done 
for an isolated cluster, we did not change the boundary 
condition and at each step of potential solving, we added 
the constant external field with ai inside the p function. 
This method might be only an approximation but as it is 
clear from the figures 1,3 and 4, the transition region from 
Newtonian to MONDian case is reproduced reasonably well 
by our method. 

Analytical solutions exist only for some special cases 
that can be subdivided as follows: 

1 - If ai ^ ao or Oe ^ ao then the value of the interpolating 
function is equal to one and the system is in the Newtonian 
regime and the velocity dispersion is given by equation ((H}. 

2 - If ae ^ ai ^ ao, the system is in the deep MOND regime 
and the external field can be neglected (isolated system). 
In this case, the line-of-sight velocity dispersion is given by 
Equation Q, which can be re- written as 



'^LOS,Ml 

GM 



= 0.049ao. 



(11) 



3 - If ai ^ <C ao, the system is externally field domi- 
nated. p(x) becomes in this case /i(ae/ao) = Oe/ao = const 
(saturation of the p function) and the system goes to a 
quasi-Newtonian regime but with an effective gravitational 
constant Ge = Gp~^ (ae/ao) — Gao/ue that is larger than 
the standard Newtonian one. Using equation ([3)l, the line- 
of-sight velocity dispersion is therefore equal to. 



= 0.035 —) ah. 

GM Oe 



(12) 



Comparing equations (|8]), (|lip and (|12p with each other sug- 
gests that the line-of-sight velocity dispersion in the general 
case should be given by 



„4 

'^LOS 



GM 

where f{ah) ~ an if ah » ao and f{ah) ~ ao if an < ao 



(13) 
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In this paper, we attempt to investigate the universal 
functional form for f{ah) for systems with a wide range of 
internal and external accelerations. 



3 N-MODY CODE 

In order to numerically solve the non-linear MOND field 
equations, recently two N-body codes hav e been developed 
ijCiotti et al. 2009 : IXiret fc Combes 20071 ). In the present 
work we apply the N-MODY code developed by the 
first group, which can be used to do numerical exper- 
iments in either MONDian or Newtonian dynamics. N- 
MODY is a parallel, three-dimensional particle-mesh code 
for the time-integration of collision-less N-body systems 
l|Londrillo fc Nipoti 200i) . The potential solver of N-MODY 
is based on a grid in spherical coordinates and is best suited 
for modeling isolated systems. N-MODY uses the leap-frog 
method to advance the particles. The code and the potential 
solver have been presented and tested by Ciotti et al. (2006) 
and Nipoti et al. (2007). 

In the present study we used a spherical grid (r, 6, ip) 
made of Nr x Ne x = 64 x 64 x 128 grid ceUs for the inte- 
gration. We use twice as many cells in the ip direction since (p 
runs from < ip < 2n while 6 runs only from < 9 < n . The 
total number of particles was in the range Np = 10^ - 10^ 
The details of the scaling of the numerical MOND models 
and code units are discussed in Nipoti, Londrillo & Ciotti 
(2007). In order to include the EFE for non- isolated systems, 
we changed the N-MODY code and put a constant exter- 
nal field in the MONDian potential solver. We also chose 
a = ai + SLe within the interpolating function as the total 
acceleration of particles. 

In the present work, the Plummer model (Plummer 
1911) was used as the initial cluster model. It has a den- 
sity distribution 



p{r) = 



3M 
4^ 



-5/2 



(14) 



where M is the total mass and rpi is the 'scale radius'. The 
half-mass radius of a Plummer model is Rh — 1.305rp; and 
the virial radius is Rv = ^rpi. The total potential energy, 
= If used in equation ([6} to calculate Vg. 



4 RESULTS 

In this Section, we present N-MODY solutions for stellar 
systems that are both isolated and non-isolated, allowing 
for different values of the external field. 



4.1 Isolated systems 

We have performed a large set of dissipationless N-MODY 
computations for isolated systems. Since the modeled sys- 
tems are in equilibrium in the Newtonian case, in the MON- 
Dian case, they initially collapse. In order to have a MON- 
Dian equilibrium initial system, we rescaled the velocities 
by an amount given by our fitting formulae ( Section 4.2 
Equations 15 ) to prevent collapse. In order to create the 
initial condition, there is another substantial method devel- 
oped by Nipoti et al. (2007b, 2008), in which the distribution 
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-1 1 

Log (a^/a„) 



Figure 1. Line-of-sight (LOS) velocity dispersion profiles for iso- 
lated stellar systems as calculated by N-MODY. In order to have 
different internal half-mass accelerations a/j, several cases with 
different half-mass radii were calculated. As expected from equa- 
tion H13| l. all curves follow the same functional form. The dashed 
line shows the asymptotic behavior in the Newtonian (equation 
[Sjl regime. The solid line shows the asymptotic behavior in deep 
MONDian (equation regime. For high internal acceleration 

(a^i 2> do): the models are consistent with the Newtonian result 
and for low acceleration (a^j <g; ao), they are consistent with the 
deep MONDian prediction. {Log = logio). 
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Figure 2. Fit of our best fitting curve (equation JTSJ) to the nu- 
merical solution for an isolated stellar system. The difference of 
each point from the fit function is presented in the inset. The av- 
erage residual of this function from our numerical solution, which 
is defined as A = \ ftheory — ffit\, is less then 10~^. {Log = logio). 
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Figure 3. Effect of different clioice of interpolation function on 
the line-of-sight velocity dispersion for systems with different in- 
ternal accelerations. Both functions have the same value in the 
Newtonian and the deep MONDian regime. In the transition zone, 
the simple function, fi2, produces a larger velocity dispersion. 
This means that if the observed velocity dispersion of a stellar 
system shows a value smaller than the MONDian prediction with 
the standard interpolation function, fii, the simple function fi2 
would not help to decrease this discrepancy. The largest difference 
between both functions is of order 20% and occurs at a^i = ag. 
{Log = logio). 

function is obtained numerically with an Eddington inver- 
sion with the far field logarithmic behavior of the MOND 
potential. Here we have used our method to set up MON- 
Dian initial condition. This method could generalis to non 
isolated systems easily (see section 4.2). 

As discussed in Nipoti et al. (2007a), all simulations 
with different masses but with the same value of ah are 
identical, in the sense that they can be simply rescaled to 
different masses, provided that M/r\ = 2ah/G remains con- 
stant. As a consequence, systems of any mass with an in the 
explored range follow the same functional form. Therefore, 
we consider only one simulation for given a^ - In order to pro- 
duce different internal acceleration regimes, we changed the 
half-mass radii of the system from 1 pc to 1 kpc. The models 
are evolved for several crossing times to reach the equilib- 
rium state, which is identified by stationary Lagrange radii 
(e.g. FigO. 

The resulting global velocity dispersions as a function 
of internal acceleration of the stellar systems are plotted 
in Fig. [T] As expected from equation (|13|l , all of them fol- 
low the same functional form. The dashed line shows the 
Newtonian prediction for the velocity dispersion (equation 
((SJ). The asymptotic behavior of the models in the Newto- 
nian regime are compatible with this analytical prediction. 
The solid line shows the analytical velocity dispersion in the 
deep MONDian regime (equation (|lip V In the low accelera- 
tion region, the numerical solutions are compatible with the 
analytical formula. At = ao, the difference between the 



Figure 4. The line-of-sight velocity dispersions for stellar sys- 
tems with an external field of = O.Olao with different internal 
half-mass-radii accelerations as calculated by N-MODY (black 
line with open squares). Ml refers to equation mil l which is for 
the isolated system and M2 refers to the quasi Newtonian case 
dominated by the external field (equation (I12l l). Due to the sat- 
uration of the ^ function in the external field dominated regime, 
the velocity dispersion curve (open squares) starts to fall in a 
quasi Newtonian way (blue line with open circles) with decreas- 
ing and deviates from the prediction of MOND for isolated 
systems. Ml, (green solid line with closed diamond). Moving to- 
wards decreasing a^, the first transition occurs near ao, 
when the system enters into the MONDian regime and the veloc- 
ity dispersion deviates from the Newtonian prediction (red dashed 
line with filled circles). The horizontal axis gives the Newtonian 
half-mass-radius acceleration. Since the Newtonian internal ac- 
celeration of a system is the square of the MONDian acceleration 
(qm = .yoArao), the point logio{ ^'" ) = —4, corresponds to 
o-h,M = O.Olao in MOND. The velocity dispersion remains on the 
horizontal line which corresponds to the isolated system until the 
internal acceleration reaches a; ae ~ O.Olao. (Log = logio). 

numerical model and the MONDian prediction is about 0.2 
in logio, which means that a{ah) — 1.3 x cjlos,mi- 

We now try to find an expression for a function /o(a;) 
where a; = — which fits the numerical results. In the New- 
tonian regime [x ^ 1), the function /o has to approach 
/o(a::) = x + const, while in the deep MONDian regime 
(a; <C 1), /o has to be constant. We therefore make an ansatz, 

/o(a::) = aln(exp(^) +6) -l-c, (15) 

for the function /q. The best-fitting coefficients are then de- 
termined by a least-squares fit to the data and are found 
to be a = 0.3314, b = 1.78, and c = -1.48. This function 
is shown in Fig. [5] as a solid line. The average residual of 
this function from our numerical results, which is defined as 
A — \ f theory — ffit\, IS less thcu 10~^. Therefore, for any iso- 
lated system, if the internal half-mass-radius acceleration, 
ah, can be measured, it is possible to find out the MON- 
Dian prediction by this function. This is especially useful 
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Log (a,/a„) 

Figure 5. External field effect on predicted line-of-sight velocity 
dispersions for stellar systems with different internal accelerations 
as calculated by N-MODY. The x-axis gives the Newtonian in- 
ternal acceleration of the system. In order to see the transition 
regime we assume several values of Og. When the internal ac- 
celeration of the system decreases, there are two transitions in 
the velocity treatment. The first transition occurs near ag 
from Newtonian into MONDian regime and the second transi- 
tion from the MONDian to quasi Newtonian regime occurs when 
the internal acceleration becomes equal to the external accelera- 
tion. The functional form of each fit curve is given in Table ((TJ- 
(Log = logio). 



for the intermediate case which has no analytical prediction 
in MOND. The corresponding formula for the velocity dis- 
persion is 

logio(az,os) = 0.25{0.331 ln[exp ( ^'^^^^^ ] + 1.78316) 

-1.48 + logio(GAfao)}. 

A simple relation exists between the three-dimensional half- 
mass radius and easier to observe two-dimensional, pro- 
jected half-mass radius Rhp : Rhp = 'yRh with 7 ~ 0.74, 
which can be used in this formulae. 

We briefly also discuss the influence of a different jj, 
function on the results. In Fig. [3] we plot the velocity dis- 
persion for an isolated system using the simple interpolation 
function fi2, and compare it with the results obtained for /ii. 
Since the simple function has a stronger MONDian effect, 
the velocity dispersion is higher than that of the standard 
function. The difference of the velocity dispersion between 
both functions at = ao is about 20%, so in order to deter- 
mine the fi function from observations, one needs to measure 
the mass and overall structure of a stellar system very ac- 
curately. As expected, in the extreme limit of <^ ao or 
3> tio, both functions predict the same value for the ve- 
locity dispersion. 



4.2 Non isolated systems 

Systems relevant for testing MOND (e.g. globular clusters or 
dwarf galaxies) usually move through the gravitational field 
of a host galaxy. Therefore, the internal dynamics is often 
influenced by the host galaxy due to the EFE of MOND. 
We assume that coriolis forces that arise in the rotating 
reference frame of the cluster and that tidal forces arising 
from a gradient of the external field can be neglected. We 
believe this to be a a good first approximation that allows 
us to focuss on the effects of the (constant) external field, 
therewith allowing us to for the first time venture into the 
intermediate MOND regime in order to study the external 
field effect numerically. 

As an example which shows the EFE on the predicted 
LOS velocity dispersion for non-isolated stellar systems with 
different internal accelerations, we choose an external accel- 
eration of fle = O.Oloo. This corresponds to a cluster or 
dwarf galaxy being at a distance of about IMpc from the 
Galactic center for an enclosed Milky Way mass of lO^^M©. 
The resulting velocity dispersion as calculated by N-MODY 
is shown in Fig. 3] The first transition occurs near ~ ao, 
when the systems enter the MONDian regime and the veloc- 
ity dispersion deviates from the Newtonian prediction. The 
velocity dispersion remains close to the MOND prediction 
for the isolated case until the internal acceleration reaches 
fli ~ ~ O.Olao at which point a second transition occurs. 
Due to the saturation of the fi function in the external field 
dominated regime, the velocity dispersion curve falls in a 
quasi Newtonian way if < fle and therefore deviates from 
the prediction of MOND for isolated systems. 

Note that the internal acceleration shown in Fig.|4]is the 
Newtonian internal acceleration, a^, = §^ , of the system. In 

the deep MOND regime, since the Newtonian acceleration 
is the square of the MONDian acceleration {um = y^ojvoo), 
the point a^^N = O.OOOlao corresponds to a^.Af ~ O.Olao. 

In order to see the effect of different external acceler- 
ations, the MONDian velocity dispersion as a function of 
internal acceleration is plotted in Fig. [S] from a weak to a 
strong external field. Note that the transition point is de- 
termined by the strength of the external field. For a smaller 
external field, the transition point occurs at a smaller accel- 
eration. As predicted by theory, for a strong external field 
(oe ^ ao) the system is completely in the Newtonian regime, 
even for a low internal acceleration. 

In order to find out the best functional form for the 
velocity dispersion as a function of the strength of the ex- 
ternal field, we use the same procedure as in the isolated 
case, but take into account the different asymptotic behav- 
ior in Fig. [5l In the Newtonian regime, {x ^ 1), f{x) still 
changes as f{x) = x + const. In the deep MONDian regime 
(x <^ 1), the systems are in the quasi-Newtonian regime and 
f{x) is again linear with the same slope as in the Newtonian 
regime. In the intermediate regime, in which the system is 
internal-acceleration dominated, af^^/GM has to be con- 
stant. A general function satisfying all these constraints is 
therefore given by 

fix) = fo{x) - aln(exp(--) + b) + c. (17) 
a 

The coefficients a, b and c depend on the external accel- 
eration. Table [1] gives their values for several values of the 
external acceleration a^. Since the asymptotic value of dif- 
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Figure 7. Upper panel: Density profile of Pal 14 for different 
cluster masses obtained by N-MODY computations. The shapes 
of the profiles are the same for all masses. Lower panel: Surface 
density profile of Pal 14 for masses as in the upper panel scaled to 
the level of data. The surface density profile shapes compare well 
with the observed density profile of Pal 14 (blue dots) as traced 
by giant stars (Hilker 2006). The meaning of the different lines is 
as in the upper panel. {Log = logio). 



ferent interpolating functions is the same, the choice of the 
/i-function does not affect systems which are in the low accel- 
eration regime. However for the higher acceleration systems 
{ah ~ do), the /i-function plays a more important role. Equa- 
tion (|17[l will for example allow it to test MOND against the 
observed global velocity dispersions of dwarf galaxies. The 
relevant external acceleration can be determined by interpo- 
lating between the points computed with N-MODY (Table 



Figure 8. Linc-of-sight velocity dispersion for Pal 14 for vari- 
ous masses as found by N-MODY. In order to compare with the 
real cluster (observational velocity dispersion of Pal 14) , the half 
mass radii of all models are fixed at 33 pc. The analytical predic- 
tions for different limiting cases are also plotted for comparison. 
Ml refers to equation (llll l. (isolated systems) and Af2 refers to 
the quasi-Newtonian case (external-field dominated case). For low 
masses, which means low internal accelerations, the prediction of 
M2 is compatible with the numerical solution. As the mass in- 
creases, the internal acceleration grows and the system enters the 
internally-dominated regime and the numerical solutions get close 
to the Ml prediction. 

4.3 The velocity dispersion of Galactic globular 
clusters 

In order to decide whether MOND or dark matter is the right 
theory to explain the dynamics of the universe, it is desir- 
able to study MOND for objects in which no dark matter is 
supposed to exist and where the characteristic acceleration 
of the stars is less than the MOND critical acceleration pa- 
rameter uq. GCs are a perfect candidate since they are the 
largest viriali zed structure that does not contain dark matter 

m. 

oore 19961 ). and their internal accelerations can be lower 
than no. Hence, G Cs may provide a good l aboratory to test 
the law of gravity (jBaumgardt et al. 20051 ). 

We choose the globular cluster Pal 14, for which there 
is a current observational effort to determine its velocity dis- 
persion (Jordi et al. 2008). We initially choose a Newtonian 
equilibrium Plummer model initially. While the half-mass 
radius of Pal 14 is about 33 pc (Hilker 2006), the mass is 
not actually known, but an observing campaign is underway 
to constrain it (Hilker et al. 2008). We change the mass in 
the wide range from [10^ — W'^]Mq and consider the half- 
mass radius to be constant. We perform numerical modeling 
to obtain the mean velocity dispersion as well as the density 
profile and velocity dispersion profile. 

Since the modeled systems are in equilibrium in the 
Newtonian case, in the MONDian case, they initially col- 
lapse and Rh is decreased before the systems virialize again. 
In order to have a MONDian equilibrium initial system, 
we increased the velocity of the Newtonian system by an 
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Figure 6. Upper left: Evolution of Lagrangian radii for a low mass cluster with stellar velocities corresponding to Newtonian virial 
equilibrium but in the deep-MOND limit. After rapid collapse, the system reaches an equilibrium state. Upper right: By increasing 
the initial velocities of particles by a factor of 2.8, calculated from equation (15), the system reaches equilibrium in MOND without 
collapsing. Lower panel: Evolution of the Lagrangian radius for a high mass cluster in the intermediate MOND regime. In the left panel, 
the velocities are not adjusted and the system still collapses. After adjusting the velocities by a factor of 1.38, calculated from equation 
(15), there is no collapse (right panel). 



amount given by our fitting formulae (Equations (15) and 
l|17p l to avoid a collapse. For low mass systems that are 
in the deep-MOND regime (a; <ti fflo) the increase is larger 
than for massive systems that are in the intermediate regime 
[tti ~ ao). In Fig. [S]we plot the evolution of the Lagrangian 
radii for two clusters with the same half-mass radius and 
different mass in the deep-MOND and intermediate MOND 
regime. In deep-MOND (low mass cluster) , after a rapid col- 
lapse, the system reaches an equilibrium state. By increasing 
the initial velocity of the particles, the collapse can be pre- 
vented. 

In Figs. [7] and |8] we show the numerical solution for 
Pal 14. We plot the density profile of Pal 14 for difi'erent 



masses and compare it with the observed profile in Fig. [7] 
(observational data from Hilker (2006)). The shape of the 
density profiles is the same for all masses, but the central 
density differs significantly. All calculated surface density 
profiles compare well with the observed density profile. It 
should be mentioned that the full observed density profile 
is not known for Pal 14 and that the surface density profile 
shown in Fig. [7] is based only on giant stars. 

Fig.[8]shows the numerical solutions for the line-of-sight 
velocity dispersion and compares it with Milgrom's analyti- 
cal predictions for the extreme limits (to see how analytical 
predictions differ from the numerical solution). Ml refers to 
equation (|ll|l which is for isolated systems and M2 refers 
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External acceleration 


a 


b 


c 


ttg = O.Olao 


0.489 


2242 


-7.694 


ae = O.OSao 


0.496 


292 


-5.680 


ae = O.lao 


0.495 


35.32 


-3.585 


ae = 0.3ao 


0.342 


8.38 


-2.119 


ae = l.Oao 


0.281 


1.00 


-0.006 


ae = lO.Oao 


0.378 


0.56 


0.477 



Table 1. Best fitting coefficients for the velocity dispersion pre- 
dicted by N-MODY simulation for various values of the exter- 
nal acceleration. The general form is given by equation (1171 , 
where fo{x) is due to an isolated cluster (equation II15I I). In case 
ae = lOoo, the function is nearly linear and the best fit function 
can also be obtained by f{x) = x — 1.467, i.e. the Newtonian 
prediction. 

to the quasi-Newtonian case which is for the external field 
dominated case (equation (|12p '). As expected, the analyti- 
cal estimates are consistent with the numerical solution ei- 
ther in the external field dominated (small mass) or internal 
field dominated case, but have significant deviations in the 
intermediate regime. For Pal 14 the external acceleration 
of the Galaxy is about ae ~ O.lao. For low cluster masses 
which means low internal accelerations, the prediction of 
M2 is compatible with the numerical solution. As the mass 
increases, the internal acceleration grows and the system 
enters the internally dominated regime and the numerical 
solution gets close to the Ml prediction. 

The velocity dispersion profiles of Pal 14 for different 
cluster masses are also plotted in Fig. [O] The velocity dis- 
persion changes slightly (about 10%) from the center to the 
half-mass-radius of the cluster. 

We would finally like to mention that the line-of-sight 
velocity dispersion in the direction of the external field is 
almost 5% lower than that perpendicular to the external 
field. Such a small difference would be unobservable. This 
can be understood as arising from the external field negating 
MONDian gravity in the direction of the field. Also, the 
cluster is elongated by a few percent along the vector Rg. 
However for larger systems such as galaxies, the anisotropy 
could be observable for a case in which ^ Oe -C oq. 



5 CONCLUSIONS 

In this work we have calculated global line-of-sight velocity 
dispersions of stellar systems in MOND for both isolated 
and non-isolated stellar systems. The velocity dispersion of 
stellar systems in MOND was so far only known in the case 
of the deep MONDian limit where all accelerations are much 
smaller than the critical acceleration, ao, and even in this 
case only if either the internal acceleration is much larger 
than the external acceleration or the internal acceleration is 
much lower than the external acceleration. We used the N- 
MODY code to calculate for the first time the hne-of-sight 
velocity dispersions of stellar systems also for the interme- 
diate regime. 

We have obtained a large set of dissipationless N- 



MODY numerical solutions for isolated systems with masses 
in the range lO'^M© to 10^ and with the Plummer model 
as the initial condition. In order to produce different inter- 
nal acceleration regimes, for each mass, we changed the half- 
mass radius of the system. We deduce the analytical formu- 
lae for the velocity dispersion of a stellar system as a func- 
tion of its half-mass-radius-internal-acceleration, ah (equa- 
tion (|13|) '). and investigate the universal functional form for 
the velocity dispersion of isolated systems (equation (I15|l '). 

We have also studied the effect of a different choice of 
the interpolation function on the line-of-sight velocity dis- 
persion for systems with different internal accelerations. We 
found that the simple function suggested by Famaey and 
Binney (2005) produces a larger velocity dispersion than 
the prediction with the standard interpolation function sug- 
gested by Milgrom (1984), with the maximum difference oc- 
curring at ah ~ ao and being of order 20%. 

Since most stellar systems (e.g. globular clusters or 
dwarf galaxies) are not isolated and usually move through 
the gravitational field of a host galaxy, the internal dynam- 
ics is often influenced by the host galaxy due to the exter- 
nal field effect (EFE) of MOND. Therefore, we have inves- 
tigated non-isolated systems, adding the external field to 
N-MODY. Our simulations reproduce previous analytic es- 
timates for stellar velocities in systems in the deep MOND 
regime (a^, ae <^ ao), where the motion of the stars is either 
dominated by internal accelerations {at 3> fle) or external 
accelerations (a^ 3> a^). In addition, we calculate the line- 
of-sight velocity dispersion for intermediate cases and derive 
for the first time analytic formulae for the line-of-sight veloc- 
ity dispersion in the intermediate regime (a; ~ a^ ~ ao) and 
found a smooth functional form for the velocity dispersion 
of stellar systems under the EFE. These formulae will allow 
to test MOND more thoroughly than was hitherto possible. 

We finally calculated the velocity dispersion of the glob- 
ular cluster Pal 14, and will compare it with observational 
data in a forthcoming paper (Jordi et al. 2008). An addi- 
tional observational study in order to constrain the mass of 
Pal 14 is also underway by our team (Hilker et al. 2008). 
In a future contribution we will also discuss the fascinating 
possibility of "freezing" a cluster on a highly eccentric or- 
bit: as a cluster moves from the Newtonian regime (small 
(jjv) to the MONDian regime on a time scale comparable 
to or faster than the internal crossing time it will retain a 
Newtonian velocity dispersion (Haghi et al. 2008). 

Additional observational efforts to determine the veloc- 
ity dispersion of stellar systems such as GCs or dSph satel- 
lites would be highly important as such data also provide 
a strict test of MOND. On the other hand, if we believe 
in MOND, these observations could be used to constrain 
the external field and consequently to put constrains on the 
potential in which the systems are embedded. Moreover it 
would be worthwhile to observe the stellar system in the 
intermediate regime to constrain the /i-function and ao. 
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Figure 9. Linc-of-sight velocity dispersion profiles of Pal 14 for different cluster masses obtained by N-MODY computation. The 
horizontal blue (dashed) line indicates the mean global velocity dispersion and the vertical red (dashed) line indicates the half-mass- 
radius of Pal 14. The velocity dispersion changes only slightly (about 10%) from the center to the half-mass-radius. 
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